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1 Introduction 
> 

CO 

^sf-N ■ The first three years of running of the LHC have probed jet physics in new ways, investigating 

previously unexplored kinematic regions. While next-to-leading-order (NLO) QCD calculations, 

00 ■ supplemented with nonperturbative corrections and parton showers, are able to describe well 

-^j- \ inclusive jet spectra over a wide range of transverse momenta extending from 20 GeV to 2 TeV, 

^^ ■ several features of LHC jet data challenge the theory. This applies in particular to the behavior 

CO I of cross sections with increasing rapidity; to correlations of multiple jets in rapidity, azimuthal 

angle, transverse energy; to non-inclusive observables probing the structure of high multiplicity 

final states. 

This article focuses on aspects of jet production which, despite the presence of a high 

transverse momentum scale, are sensitive to soft gluon processes and QCD infrared physics. 

C^ ■ We start with inclusive cross sections in Sec. [2] and discuss the role of parton showering and 

nonperturbative effects in the context of matched NLO-shower event generators. In Sec. [3] 

we consider forward jets, and examples of multi-jet correlations. Sec. 2] examines 6-flavor jets. 

Sec.[5]takes a further look at jet correlations from the viewpoint of multiple parton interactions, 

emphasizing the role of energy flow variables. Sec. [6] addresses motivation and prospects for 

extending jet measurements to lower transverse momenta than is presently done. 

2 Inclusive jet production 

Measurements of inclusive jet production are carried out at the LHC [1, 2j over a kinematic 
range in transverse momentum and rapidity much larger than in any previous collider exper- 
iment [3]. Baseline comparisons with Standard Model theoretical predictions rely either on 
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next-to-leading-order (NLO) QCD calculations, supplemented with nonperturbative (NP) cor- 
rections [HE] estimated from Monte Carlo event generators, or on NLO-matched parton shower 
event generators |4l[5]. The upper panels in Fig.[T][l report the first kind of comparison, show- 
ing that the NLO calculation agrees with data at central rapidities, while increasing deviations 
are seen with increasing rapidity at large transverse momentum pt P. The question arises 
of whether such behavior is associated with higher-order perturbative contributions or with 
nonperturbative components of the cross section. 
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Figure 1: Inclusive jet spectra fJ]/ compared with (top) NLO+NP results and (bottom) NLO- 
matched shower results. 



The lower panels in Fig. [T] [1^ show the second kind of comparison based on Powheg 
calculations [6 , in which NLO matrix elements are matched with parton showers [IlIH]. This 
improves the description of data, indicating that higher-order radiative contributions taken 



into account via parton showers are numerically important. At the same time, the results show 
large differences between Powheg calculations interfaced with different shower generators, 
Herwig [7^ and Pythia [8j, in the forward rapidity region. This region is sensitive to the 
details of parton showering corrections. 

These observations make it apparent that QCD contributions beyond NLO are essential 
for the understanding of LHC jet data. In the following subsection we discuss showering and 
nonperturbative effects. 

2.1 Nonperturbative and showering corrections 

Using leading-order Monte Carlo (LO-MC) generators [3 [8], the nonperturbative correction 
factors are schematically obtained in (TJ [2] as 

^0 —^^ LO-MC /^^ LO-MC ' \^) 

where {ps + mpi -\- had) and {ps) mean respectively a simulation including parton showers, 
multiparton interactions and hadronization, and a simulation including only parton showers in 
addition to the LO hard process. 

While this is a natural way to estimate NP corrections from LO+PS event generators, it is 
noted in [9] that when these corrections are combined with NLO parton-level results a potential 
inconsistency arises because the radiative correction from the first gluon emission is treated at 
different levels of accuracy in the two parts of the calculation. To avoid this, Ref. [9] proposes 
a method which uses NLO Monte Carlo (NLO-MC) generators to determine the correction. In 
this case one can consistently assign correction factors to be applied to NLO calculations. This 
method allows one to study separately correction factors to the fixed-order calculation due to 
parton showering effects. To do this, Ref. [9 introduces the correction factors K^^ and K^^ 
as 

T^NP _ j^Aps+mpi+had) , i\Aps) /r)\ 

^ — ^^ NLO-MC / ^^ NLO-MC ' V^/ 

^ — ^^NLO-MCl ^^ NLO-MC ' W/ 

where the denominator in Eq. (|3]) is defined by switching off all components beyond NLO in 
the Monte Carlo simulation. 

The factor K^^ in Eq. (|2]) differs from K^^ because of the different definition of the hard 
process. In particular the multi-parton interaction p^ cut-off scale is different in the LO and 
NLO cases. Numerical results are shown in Fig. [2l The factor K^^ in Eq. ([3]), on the other 
hand, is new. It singles out contributions due to parton showering and has not been considered 
before. Unlike the NP correction, it gives finite effects also at large pt- The results in Fig. [3] 
show in particular that this correction is not merely a rescaling but is y and pt dependent, 
especially when rapidity is non-central. 

The result in Fig. [3] comes from initial-state and final-state showers. These are interleaved 
so that the combined effect is nontrivial and cannot be obtained by simply adding the two |9]. 
In general the effect from parton shower is largest at large |^|, where the initial-state parton 
shower is mainly contributing at low j9t, while the final-state parton shower is contributing 
significantly over the whole pt range. It is observed in [9^ that the main initial- state showering 
effect comes from kinematical shifts in longitudinal momentum distributions [10], due to com- 
bining collinearity approximations with the Monte Carlo implementation of energy-momentum 
conservation constraints. The effect of the kinematic shifts is illustrated in Fig. 2] |9j, showing 
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Figure 2: The NP correction factors to jet transverse momentum distributions obtained from 
Eq. (Cp and Eq. (0), using Pythia and Powheg respectively, for \y\ < 0.5 and 2 < \y\ < 2.5. 
Left: R = 0.5; Right: R = 0.7. JEj 
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Figure 3: The parton shower correction factor to jet transverse momentum distributions, ob- 
tained from Eq. (0) using Powheg for \y\ < 0.5 and 2 < \y\ < 2.5. Left: R = 0.5; Right: 

R = o.7. m 



the distribution in the parton longitudinal momentum fraction x before parton showering and 
after parton showering. We see that the longitudinal shift is negligible for central rapidities but 
becomes significant for y > 1.5. 

In summary, the nonperturbative correction factor K^^ introduced from NLO-MC in 
Eq. ([2]) gives non-negligible differences compared to the LO-MC contribution [U [2] at low 
to intermediate jet pt, while the showering correction factor K^^ of Eq. (|3]) gives significant 
effects over the whole pr range and is largest at large jet rapidities y. Because of this y and 
Pt dependence, taking properly into account NP and showering correction factors changes the 
shape of jet distributions, and may thus influence the comparison of theory predictions with 
experimental data. We anticipate in particular that taking account of the showering correction 
factor will be relevant in fits for parton distribution functions using inclusive jet data. 
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Figure 4: Distributions in the parton longitudinal momentum fraction x before (POWHEG) and 
after parton showering (POWHEG-hPS), for inclusive jet production at different rapidities for 
jets with Pt > 18 GeV obtained by the anti-kt jet algorithm /77]/ with R = 0.5. Shown is the 
effect of intrinsic kt, initial (IPS) and initial -h final state (IFPS) parton shower. 0/ 



3 Forward jets 



A program of jet physics in the forward region can be carried out for the first time at the 
LHC [12], by exploiting the large phase space opening up at high center-of-mass energies, and 
the unprecedented reach in rapidity of the experimental instrumentation [13]. First forward jet 
measurements have been performed by LHC experiments [14^ J5l [16] . 



The evaluation of QCD theoretical predictions for forward jets is made complex by high- 
er production occurring in a region characterized by multiple hard scales, possibly widely 
disparate from each other [12] . In addition to the parton-shower kinematic effects discussed in 
Sec. 2, when jets are measured at large separations in rapidity dynamical contributions from 
soft multi-gluon radiation ^7\ [181 [E] set in, calling for perturbative resummations of large- 
rapidity logarithms to all orders in a^. This motivates current studies based on the BFKL 
equation [20l [211 [221 [23]. Furthermore, with increasing rapidities the nonperturbative parton 
distributions are probed in highly asymmetric kinematic regions near the boundaries x ^ and 
X ^ 1 of partonic phase space [121 [IS] • This in turn implies that contributions from multiple 
parton collisions [24l [25] are potentially enhanced [26l [271 [21] • Effects of multiple scatterings 
may be studied in di-jet as well as di-hadron spectra [29], in pp collisions and in collisions of 
hght [30] and heavy [SB [32l [33] nuclei. 

Forward jets enter the LHC physics program in an essential way both for new particle 
discovery processes (e.g., Higgs searches [34 in vector boson fusion channels, jet studies in 
decays of highly boosted heavy states [35 ) and for new aspects of standard model physics 
(e.g., small-x QCD and its interplay with cosmic ray physics [36l [37], studies of high-density 
parton matter [38l |39]). The first forward jet measurements [HI [151 [IS] show that, while 
inclusive forward jet spectra are roughly in agreement with predictions from different Monte 
Carlo simulations, detailed aspects of production rates and correlations are not well understood 
yet. Recent phenomenological analyses are carried out in [4Ql[4T]. 
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Figure 5: Ratio theory /data ml for dijet events with a central and a forward jet as a function 
of the forward jet transverse momentum. 



An example is given by the di-jet observables [42 associated with events containing a for- 
ward and a central jet. Experimental measurements and Monte Carlo comparisons are shown 
in Fig. [5] [14]. The results indicate that none of the Monte Carlo generators describes the 
data well in all regions; in particular NLO-matched calculations from Powheg give large dif- 
ferences in the forward jet pt distribution when combined with different parton showers, see 
Powheg+Herwig vs. Powheg+Pythia. 

In [321 [l3j this behavior is investigated by studying Ai? jet distributions designed to provide 



a measure, in azimuth and rapidity space, of the extent to which jets are dominated by hard 
partons in the matrix element or originate from showering. Large contribution to jets from 
showering are found [42] in the case of asymmetric parton kinematics, i.e, when one of the initial- 
state showers goes down to small x. Ref. [l3] furthers this study by considering the central jet 
transverse energy spectrum, in di-jet events with a central and a forward jet, using the NLO 
event generator Powheg matched with parton showers Pythia and Herwig. Fig. [6] [43 shows 
results for the two cases, normalized to the result obtained by switching off parton showering. 
The marked differences between the two cases are consistent with the findings in [14 , and with 
the large contribution to jets from showering found in [12]. In particular in the forward-central 
events considered high-rapidity correlations turn out to affect the behavior of jet distributions 
in the central region. 
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Figure 6: Ratio of NLO-hshower to no- shower results for dijet events with a central and a 
forward jet as a function of the central jet transverse momentum {43^ . (left) Et > 10 GeV; 
(right) Et >30 GeV. 



A classic test of QCD high-energy resummation for jets at large rapidity separations [17] 
is given by the azimuthal decorrelation between jets. Fig. [7] [12] shows the cross section as a 
function of the azimuthal distance Acj) between central and forward jets reconstructed with the 
Siscone algorithm [44] {R = 0.4), for different rapidity separations. It shows results computed 
by Pythia Monte Carlo [IS], with and without multi-parton interactions, and by Cascade 
Monte Carlo [46^, which includes small-x gluon coherence effects |47J in the initial-state shower. 
The main point is that the decorrelation as a function of Arj increases in Cascade as well as 
in Pythia, respectively as a result of finite-angle gluon radiation in single-chain parton shower 
or as a result of multiple-chain collinear showers; but while in the low Et region (Fig. [7] (left)) 
this is similar between Cascade and Pythia with multiparton interactions for A77 < 4, in 
the higher Et region (Fig. [3 (right)) the influence of multiparton interactions in Pythia is 
small and Cascade predicts everywhere a larger decorrelation. We will come back to this and 
discuss correlations further in Sec. 5 from the point of view of energy flow observables. 

While the specific results shown in this section refer to forward-central jet correlations, it 
is interesting to consider extensions to the forward-backward kinematics. This will allow one 
to address the large- A?/ di-jet data sets [13 US], currently rather poorly understood; to search 




Figure 7: Cross section versus azimuthal distance Acj) between central and ^ 

ent rapidity separations Ar], for jets with transverse energy Et > 10 GeV (left) 

GeV (right) m. 




for Mueller-Navelet effects [17] ; to analyze backgrounds in Higgs boson studies i34J from vector 
boson fusion channels. In particular, one may be able to extract information on Higgs properties 
and couplings from jet kinematics [48]. In this case too finite-angle radiative contributions to 
single-chain showers, extending across the whole rapidity range, affect the underlying jet activity 
accompanying the Higgs ^ and may give competing effects to multiple-parton interactions. 



4 B-jets 

Some of the features observed for inclusive jets in Sec. 2 are also present in 6-flavor jets [50l 
[51]. Fig. [8] ^^ shows a comparison of the measured B-]ei transverse momentum spectra 
with matched NLO-shower calculations using MC@NLO [52 . The description of the data is 
generally good at central rapidities, while at large y and large pr the Monte Carlo is above the 
data. Similar behavior is shown by comparisons with Powheg [53] in [5T] . 

Fig. [9] [9] studies kinematic corrections to 5-jet production due to longitudinal momentum 
shifts in the initial-state parton shower, similar to those discussed for inclusive jets in Fig. H) 
For 5-jets in different rapidity regions |^, the gluon x distribution is plotted from Powheg 
before parton showering and after including various components of the parton shower generator, 
using the Pythia parton shower Z2 [54] (including hadronization components to identify B- 
jets). Fig. [9] shows similar shifts in longitudinal momentum with increasing rapidity as in the 
inclusive jet case. A better understanding of B production in this region will affect studies of 
the Higgs to hh decay channel, e.g. in the associated production with vector bosons. 

Although the explicit calculations in Fig. [9] are performed using a particular NLO-shower 
matching scheme (Powheg), the effect is common to any calculation matching NLO with 
collinear showers. As discussed in [9l |T0], the kinematic shifts due to the momentum reshuf- 
fling can affect predictions of matched NLO-shower calculations both through the perturbative 
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weight for each event and through the evaluation of the parton distribution functions. In calcu- 
lations using integrated parton density functions this implies that correction factors as discussed 
in Sec. 2 have to be applied after the evaluation of the cross section. On the other hand, we note 
that this is avoided in approaches using transverse momentum dependent pdfs [55l [56l [571 EH] 
from the beginning (TMDs or uPDFs). It will thus be of interest to study it quantitatively in 
Monte Carlo generators which implement these pdfs [l6l \59\ [60]. 

We note further that helicity amplitudes techniques are being developed [61] for multi-leg 
processes with up to two off-shell gluons in the high-energy limit ^8\ , and applied to complex 
final states containing heavy quarks [62], including bb + jets and bb + vector bosons. Once 
combined with TMD pdfs and parton showers, these multi-parton amplitudes can provide 
a computational tool to overcome the limitations of collinear approximations, and take into 
account both kinematical and dynamical effects associated with multiple-scale processes. 



5 Multi-parton interactions and energy flow variables 

Besides jet cross sections, event shape variables are studied at the LHC [63l [6l] and used 
to characterize the events' energy flow and the structure of multi-jet flnal states. Baseline 
descriptions of multi-jet events are obtained by methods merging parton showers and multi- 
leg hard matrix elements [5^ . First measurements of LHC hadronic event shapes [53] point to 
parton showering effects dominating over contributions from hard matrix elements evaluated 
at high multiplicity. Jet shape variables describing the jet's internal structure and the energy 
flow within a jet are also studied [65]. These observables are used in searches for potential 
new physics signals from decays of massive states in the boosted regime [35 , or away from 
it [66]. Besides jet substructure, such observables are sensitive to soft physics effects, including 
underlying events, pile- up, multiple parton interactions [26] . 
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Figure 10: Transverse energy flow fTUil in the inter- jet region: (left) particle flow; (right) mini- 
jet flow. 



Energy flow measurements [67 in minimum bias and dijet events, designed to investigate 
properties of the soft underlying event, emphasize the difficulty [68] in achieving a unifled 



10 



underlying event description from central to forward rapidities, based on Pythia [45] Monte 
Carlo tuning. Forward-backward correlations ^^ in minimum bias may help analyze the event 
structure. Complementary to the above measurements are transverse energy flow observables 
associated with the production of jets widely separated in rapidity ^, sensitive to harder 
color radiation, and useful for studies of showering and of multi-parton interactions [71J. The 
transverse energy flow may be defined by summing the energies over all particles in the final 
states above a minimum Et^ or alternatively [70^ by first clustering particles into jets by means 
of a jet algorithm, and then constructing the associated energy flow from jets with transverse 
energy above a given lower bound go- In the latter case one measures a (mini)jet energy flow, 
and infrared safety is ensured by the use of the clustering algorithm. 

Figs. [TOl and [TTI report results for the particle and minijet energy flow associated with pro- 
duction of central and forward jets [70^ from three Monte Carlo event generators: the k^-shower 
Cascade generator [46], to evaluate contributions of high-energy logarithmic corrections; the 
NLO matched Powheg generator [6], to evaluate the effects of NLO corrections to matrix 
elements; Pythia Monte Carlo [IS], used in two different modes: with the LHC tune Zl [51] 
(PYTHlA-mpi) to evaluate contributions of multi-parton interactions, and without any multi- 
parton interactions (PYTHlA-nompi). 
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Figure 11: Azimuthal dependence of the mini- jet energy flow f7U\l for different rapidity ranges: 
(left) central-jet; (middle) intermediate; (right) forward-jet. 

Fig. [10] shows the pseudorapidity dependence of the transverse energy flow in the region 
between the central and forward jets. The particle energy flow plot on the left in Fig. [TOl shows 
the jet profile picture, and indicates enhancements of the energy flow in the inter-jet region 
with respect to the PYTHlA-nompi result from higher order emissions in CASCADE and from 
multiple parton collisions in PYTHlA-mpi. On the other hand, there is little effect from the 
next-to-leading hard correction in Powheg with respect to PYTHlA-nompi. The energy flow 
is dominated by multiple-radiation, parton-shower effects. The mini-jet energy flow plot on 
the right in Fig. [TOl indicates similar effects, with reduced sensitivity to infrared radiation. As 
the mini-jet flow definition suppresses the contribution of soft radiation, the Cascade and 
PYTHIA-mpi results become more similar in the inter-jet region. Distinctive effects are also 
found in [70J by computations in the region away from the jets. 

Fig. [11] illustrates the azimuthal dependence of the mini-jet transverse energy flow. Here 
A^ is measured with respect to the central jet. The A^ distribution is shown for three different 
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rapidity ranges, corresponding to the central-jet, forward-jet, and intermediate rapidities. As we 
go toward forward rapidity, the Cascade and PYTHlA-mpi calculations give a more pronounced 
flattening of the A(/) distribution compared to Powheg and PvTHlA-nompi, corresponding to 
increased decorrelation between the jets. 

The above numerical results indicate that soft, finite- angle multi-gluon emission over large 
rapidity intervals gives sizeable contribution to the inter-jet energy flow. Non-collinear cor- 
rections to single-chain showering may in particular influence the rates for multi-parton inter- 
actions [26]. This also underlines the relevance of approaches which aim at a more complete 
description of initial state dynamics by generalizing the notion of parton distributions [56l [57] , 
both for quark-dominated [72] and gluon-dominated [73] processes. It will be of interest to 
extend energy flow studies using different triggers, e.g. vector boson + jet [74] at large rapidity, 
which may help constrain multiple interaction rates [45l [54l [75] . 

6 Towards low pt 

Hadron jets are measured at the LHC over a wide range in transverse momenta from 20 GeV to a 
few TeV. It is of interest to explore to what extent the picture of jets provided by the algorithms 
used for jet reconstruction (and corresponding shower Monte Carlo event simulations) can be 
pushed towards the semi-soft region of lower transverse momenta. As discussed earlier, low- 
Pt jets may be used to investigate QCD physics at small x. A jet-like description of low-pT 
hadronic final states may provide insight into the region where the LHC has produced striking 
evidence ^J of long-range di-hadron correlations in high multiplicity events. 

It has been observed in [77 that if jet measurements at the LHC are extended down to 
transverse momenta of the order of a few GeV one can define a visible leading jet cross section 
sensitive to the unitarity bound set by the inelastic proton-proton rate which has recently been 
measured [78l [79l [80] . This can be done within the range of acceptance of the measurement 
without using any extrapolation [77 . Because of the low transverse momenta, this will rely pri- 
marily on jets constructed from charged tracks. Given the decay in particle tracking capabilities 
with increasing rapidity, one may focus on the central pseudorapidity range. 

The main interest of these measurements is the possibility to investigate the leading jet 
cross section near the pr region, pt = O (a few GeV), where the inelastic pp production rate 
is saturated [77] • At the LHC, for the first time in collider experiments, this occurs in the 
weakly coupled region. Even though at weak coupling, dynamical effects slowing down the 
rise of the cross section in this region involve strong fields and non-perturbative physics. In 
order to analyze the saturation region, Ref. [771 introduces an event cross section, defined as 
an integral over the differential leading jet cross section, which does not depend on the jet 
multiplicity. Fig. [12] shows the visible jet cross section using Pythia with jets reconstructed 
by the anti-k^ algorithm [11] for i? = 0.5 down to low transverse momenta. 

In Fig. [12] (left) the perturbative result reaches the inelastic bound [78^ for minimum p^ — 4 
GeV. In the region just above this value, pr = 0(10) GeV, effects responsible for the taming of 
the cross section set in. Fig. [12] (right) shows the cross section based on the model [25l[8T], in 
which the rise of the cross section is tamed at small values of pr by introducing a pro cut-off 
parameter obtained from fits to describe measurements of the underlying event. 

Fig. [13] [77 shows a comparison of the jet cross section for ptq = 0, including parton shower 
and hadronization, with the cross section obtained from Pythia including the pto niodel. In 
Fig. [13] (right) we show the effect of multi-parton interactions. Especially in the region of 
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Figure 13: The cross section J771 as a function of pTmin cls predicted by Pythia in the range 
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and hadronisation, while the dashed (red) line shows the prediction with pro 7^ 0; (left) is without 
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Pt < 10 GeV a clear deviation from the pro = prediction is visible. Measuring the jet cross 
section in this range would probe the transition from the large-pT perturbative behavior to the 
medium to low pr region where (weak-coupling) unitarity corrections set in. We note that the 
analysis [82J already illustrates the feasibility of measuring jets at low pr- However, event cross 
sections such as that described above to study unitarity have not yet been examined. Results 
from different Monte Carlo models in [82] are effectively normalized to the lowest pt bin. 

Besides pp collisions, low pt jet measurements are relevant in pA and A A nuclear collisions. 
If the inelastic cross section is measured in AA and pA, they may be useful to characterize 
properties of final states in terms of jets or flows [83], and investigate transverse momentum 
dependent dynamics and multi-parton interactions [84] [85] . 

Compared to the collinearly-factorized cut-off model [25l [81] , a physically different picture 
of the turn-over region is provided by approaches based on TMD initial-state distributions [56l 
[58l [86] incorporating small-x gluon coherence ^. In this case the rise of the cross section is 
slowed down by finite transverse momenta in the initial state. The effect of combining this with 
nonlinear terms in the evolution of parton cascades [87] and the dependence of the turn-over 
region on the infrared behavior of the strong coupling [88] warrant investigations. Extending 
jet definition and measurements to the region of low transverse momenta may give us radically 
new information on the underlying physical picture of high-energy multi-particle production. 

7 Conclusion 

Jets enter in essentially all areas of the LHC physics program, from new particle discovery 
processes to precision studies of the Standard Model to searches for physics beyond the Standard 
Model. Experimental results in the first three years of running of the LHC have provided us 
not only with solid confirmation of our understanding of jet physics from QCD but also with 
challenges and surprises. These involve both the high-pT region, featuring several new effects 
of QCD parton showers, and the low-pT region, probing nonperturbative components of jet 
cross sections. Much is to be learnt about the theory of hadronic jets and its appHcations to 
precision phenomenology as forthcoming LHC data analyses and measurements keep exploring 
jets in new regions of phase space and probing different sectors of QCD. 
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